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ABSTRACT 


This  work  consists  of  the  development  of  the  six  degree 
of  freedom  non  linear  model  of  a  sea  launched  generic  bank-to- 
turn  cruise  missile  attacking  a  medium  sized  combatant  ship. 
Two  guidance  and  control  schemes  are  compared  in  the  terminal 
phase.  The  first,  or  baseline  guidance  scheme  (pop  out 
maneuver),  uses  a  50-foot  altitude  hold  for  an  ingress  phase, 
followed  by  a  pop  out  maneuver,  and  then  an  attack  phase  which 
uses  proportional  navigation  in  elevation  and  azimuth  planes 
along  with  bank-to-turn  maneuvering.  The  second  scheme 
(sea  skimmer)  uses  identical  ingress  and  attack  phases  but 
eliminates  the  pop  out  maneuver.  Miss  distances  for  both 
schemes  are  compared  while  varying  missile  roil  rate  limit, 

SCM  blinking  frequency,  and  burn  through  ranges. 
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I.  INTRODUCTION 


The  relative  merits  of  bank-to-turn  versus  skid-to-turn 
missiles  have  been  argued  for  years.  Bank-to-turn  missiles 
which  must  roll  to  a  commanded  bank  angle  before  a  lateral 
acceleration  can  be  commanded  in  the  appropriate  direction  are 
inherently  slower  to  respond  to  target  maneuvering  than 
skin-to-turn  missiles  of  comparable  lateral  acceleration 
capabilities.  This  is  particularly  true  if  the  bank-to-turn 
missile  has  limited  roll  authority  or  limited  roll  rate 
capability.  On  the  other  hand,  as  discussed  in  Gonzalez, 

{Ref.  1},  it  is  sometimes  difficult  to  build  a  skid-to-turn 
missile  which  meets  certain  operational  goals  and  performance 
criteria  and  also  has  sufficient  control  authority  in  both 
the  y  and  z  directions  to  generate  the  large  lateral  accel¬ 
erations  required  to  perform  skid-to-turn. 

In  this  thesis  the  terminal  guidance  accuracy  of  a 
bank-to-turn  cruise  missile  is  studied  for  a  sea  launched, 
sea  target  cruise  missile  of  conventional  wing-tail  config¬ 
uration.  It  is  assumed  that  roll  authority  is  extremely 
limited  due  to  operational  design  restrictions  which  impose 
the  use  of  low  authority  ailerons  or  differential  stabilator 
(no  ailerons)  for  rolling.  Lateral  acceleration  is  also 
assumed  to  be  extremely  limited  in  the  y  direction  due  to 
the  presence  of  only  a  small  vertical  surface  at  the  tail 
(vertical  stabilizer  and  rudder). 
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Two  bank-to-turn  guidance  schemes  are  compared  for 
accuracy  against  a  medium  sized  combatant  ship  employing 
electronic  countermeasures  (an  ECM  blinker  mounted  aft  on 
the  ship).  The  first  baseline  scheme  employs  a  pop-out 
maneuver  consisting  of  a  low  altitude  run  in,  azimuth  offset, 
pop-up  and  roll  into  the  target  using  proportional  navigation 
in  azimuth  and  elevation.  The  second,  or  sea-skimmer  scheme, 
employs  a  straight-in  low  altitude  attack.  Miss  distances 
for  both  schemes  are  compared  while  varying  missile  roll 
rate  limit,  ECM  blinking  frequency  and  burn  through  ranges. 

A  CSMP  III  simulation  was  coded  to  perform  the  study. 
Missile  dynamics  were  represented  by  6  degree  of  freedom 
nonlinear  equations  of  motion  with  table  look-up  aero¬ 
dynamic  coefficients.  These  coefficients  are  representative 
of  a  generic  cruise  missile  with  a  conventional  wing- tail 
configuration  of  limited  roll  authority.  Inner  loop 
augmentation  and  autopilot  modes  were  designed  to  improve 
missile  damping  in  all  axes,  provide  commanded  load  factor, 
commanded  bank  angle,  and  turn  coordination.  Autopilot 
outer  loops  were  designed  to  provide  altitude  hold,  and 
vertical  flight  path  angle  hold.  Guidance  loops  were  designed 
to  provide  proportional  navigation  in  elevation  and  azimuth. 

Since  this  study  involves  the  influence  of  flight 
dynamic  parameters  (roll  rate  limits)  on  terminal  guidance 
accuracy,  no  extensive  tracker  modeling  is  employed  in  the 
simulation.  It  is  assumed  that  the  seeker  always  tracks 


the  point  target  perfectly.  The  point  target  is,  however, 
moved  by  an  ECM  blinker  and  contaminated  with  a  glint  model. 

In  Chapter  II  the  equations  used  to  represent  missile 
dynamics  are  presented  along  with  the  methodology  used  to 
build  up  aerodynamic  forces  and  moments.  In  Chapter  III 
the  development,  design  and  simulation  of  the  missile  auto¬ 
pilot  modes  are  presented.  In  Chapter  IV  the  design  of  the 
guidance  systems  and  a  complete  description  of  the  two 
guidance  schemes  is  presented.  In  Chapter  V  the  CSMP 
simulation  is  presented  followed  by  the  results  and  conclusions 
of  this  study. 
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II.  DEVELOPMENT  AND  SIMULATION  OF  MISSILE  DYNAMICS 


In  this  chapter  the  linear  and  nonlinear  mathematical 
models  are  developed  which  are  used  to  describe  the  flight 
dynamics  of  a  generic  bank-to-turn  cruise  missile.  The 
linearized  mathematical  model  of  missile  dynamics  is  used 
to  design  the  missile  autopilot  control  laws  and  guidance 
laws.  The  nonlinear  mathematical  model  is  used  to  accurately 
represent  missile  dynamics  in  a  CSMP  computer  simulation  of 
a  sea  launched  generic  cruise  missile  attacking  a  medium 
sized  combatant  ship.  Since  the  purpose  of  this  investigation 
is  to  perform  a  detailed  evaluation  of  terminal  control 
laws,  it  was  decided  that  a  linear  dynamic  simulation  would 
not  adequately  represent  missile  motion;  hence  a  full  six 
degree  of  freedom  (6  DOF)  nonlinear  dynamic  simulation  was 
encoded.  The  overall  system  that  is  modeled  in  this  study 
is  depicted  in  block  diagram  format  in  Figure  2-1  and  will 
be  developed  in  this  thesis. 

A.  MISSILE  EQUATIONS  OF  MOTION 

The  full  nonlinear  six  degree  of  freedom  rigid  body, 
dynamic  equations  of  motion  are  used  to  represent  the  motion 
of  a  generic  bank-to-turn  cruise  missile.  The  aerodynamic 
forces  and  moments  are  built  up  from  representative  tabular 
coefficient  data  for  the  generic  missile.  These  tabular 


data  are  given  in  Appendix  A  in  graphical  form  along  with  the 
generic  missile's  physical  and  geometric  characteristics. 

1 .  Assumptions 

The  assumptions  used  in  the  development  of  the 
equations  of  motion  are  given  below. 

a.  The  earth  is  flat,  does  not  rotate,  and  is  fixed 
in  inertial  space. 

b.  The  mass  of  the  missile  is  constant. 

c.  The  missile  is  a  rigid  body. 

d.  The  mass  distribution  of  the  missile  is  constant. 

e.  Engine  angular  momentum  is  neglected. 

f.  The  missile  is  symmetric  about  the  body's  xy 
plane.  Therefore,  the  products  of  inertia  Ix^  and  I  are 
zero. 


g- 

body  x  axis 
zero. 


he  engine  thrust  line  is  parallel  to 
Thus,  the  thrust  components  T^  and  T 


the  missile 

are 

z 


h.  The  density  of  the  atmosphere  is  constant. 

i.  The  engine  thrust  line  passes  through  the  missile 
center  of  mass.  Therefore,  the  moments  due  to  thrust  L^,, 

M,j,  and  N^,  are  zero. 

j .  The  wind  is  calm. 

k.  Right  aileron  trailing  edge  down  is  positive 
deflection.  Negative  elevator  deflection  yields  nose  up 
pitching  moment  (positive). 
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2.  Coordinate  System 

An  earth  fixed  coordinate  system  is  established_with 

its  origin  fixed  at  the  initial  position  of  the  missile.  The 

X„  axis  points  to  true  north,  the  axis  points  east,  and 

the  Z-.  axis  points  toward  the  center  of  the  earth.  Altitude 
c< 

(H)  therefore  equals  the  negative  of  Z ...  The  system  is 
considered  to  be  an  inertial  system. 

A  stability  axis  coordinate  system  is  used  with 
coordinates  x,  y,  2  fixed  at  the  missile  center  of  mass  and 
oriented  such  that  the  x  axis  lies  along  the  vehicle's 
forward  velocity  vector  in  steady  state  balanced  cruise 
flight  at  0.75  Mach  number.  The  y  axis  is  chosen  perpendicular 
to  the  plane  of  symmetry  and  is  oriented  out  the  right  wing. 

The  z  axis  is  chosen  perpendicular  to  the  x  axis  in  the  down 
direction  and  in  the  plane  of  symmetry. 

3.  The  Equations  of  Motion 

The  following  nonlinear  equations  as  developed  by 
Hewett  {Ref.  2}  describe  the  motion  of  the  missile  in 
6  degrees  of  freedom. 


m  (fr-VR+WQ)  =  -mg  sin  9  +  X  +  T 

(la) 

m  (V+UR-WP)  =  mg  sin  cos  9  +  Y 

(1b) 

m  (W-UQ+VP)  *  mg  cos  cos  9  +  Z 

(1c) 

PIxx-(&tP«)lx2+  RQU^-Iyy)  =  LA 

(2a) 

«1iry+PR(Ixx-Iz,)  +  (P^-R^)IXy  -  Ma 

(2b) 

ftl„-flx,+p<»(VIxx,+,SRIx.  *  »A 

(2c) 
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<$  =  P+ (Q  sin  $  +  R  cos  $ )  tan  0 
9  =  Q  sin  $  -  R  sin  $ 

4*  =  (Q  sin  $  +  R  cos  4> )  sec  0 


.(3b) 


Equations  (la),  (1b)  and  (ic)  describe  the  transla¬ 
tional  motion  of  the  missile.  Equations  (2a),  (2b),  and 
(2c)  represent  the  rotational  motion.  Equations  (3a), 

(3b)  and  (3c)  are  the  Euler  relations  for  bank,  pitch  and 
yaw  angles. 

Equations  (4a),  (4b),  and  (4c)  describe  the  missile's 

position  referenced  to  the  earth  fixed  system  (Xg,  Yv,  ZT). 

$■ 

t-,  =  U  cos  fcos  0+  V(  cosY  sin9  sin$>  -  sin'P  cos<!> )  (4a) 

Zi 

+  W(  cosfsin©  cos$  + sinfsin*!’ ) 

£„  =  U  sin  fcos  0+  V(  sinfsin©  sin<t  +  cos4,cos$  )  (4b) 

sh 

+  W(sin4'sin0cos$-cos4'sin$ ) 

=  Usin0-Vcos0  3in$-Wcos0cos<*>  (4c) 

Ci 

Equations  (5),  (6)  and  (7)  are  also  required. 
Equations  (5)  and  (6)  yield  angle  of  attack  and  side  slip 
angle  respectively  in  terms  of  velocity  components  U,  V  and 
W.  Equation  (7)  represents  the  total  velocity  of  the  missile, 


<*  =  arctan 


=  arcsin 


W 

(U2+W2)^ 


B  =  arctan  — 3—  ^ 


=  arcsin 


(U^+V2+W2)^ 


(U2+V2+W2)^ 


Figure  2-2  illustrates  the  positive  directions  and 
locations  of  forces,  moments  and  velocities  in  the  body 
coordinate  system  (x,  y,  z). 


Aarodyry 


l 


lc  and  Thrust  Moments 


Linear  and 
Velocities 


Rotational 


(Angular] 


Figure  2-2  Body  Coordinate  System  Description 


Equations  (la),  (1b),  (1c)  and  (2a),  (2b),  (2c)  can  be 
rewritten  in  state  format  as  follows: 


ft  =  -gsine  +  VR-WQ  +  -J-  +  jjj- 

V  =  gsin$ cosQ -UR+WP  +  ■— 

W  =  gcos$cosQ+UQ-VP  +  ■— * 


*  =[la 


(I  I  -I  )  (9a) 

xx  zz  xz 


9  '  1, 


R  = 

> 


LAIx,+1,*Ixx-p,>(Iyyrxx-I«-I«)-i,,5(Ixx+IzS-Iyy)Ixz  1 

a 


yy  xx  xx  xz 


'xx  zz  yy  xz 


(I  I  -I  ) 

XX  zz  xz 


Now  integrating  equations  (3).  (4),  (8)  and  (9)  results 
in  U,  V,  W,  ?,  ;,  R,  i>,  3,  9,  X£,  Y_  and  Z£  (-H).  The 
missile's  velocities,  angular  rates,  and  angles  and  positions 
have  thus  been  described. 

4.  Trim  Equations 

For  straight  and  level  flight  at  constant  velocity 

•  ••••• 

the  following  variables  are  all  zero:  U,  V,  W,  P,  Q,  R, 


l,  5,  2S,  pgs,  Qss,  Rss,  tss»  lass’  nass'  YSS’  ySS’ 


Thus  the  trim  equations  reduce  to: 
SS 

XSStTx  =  “8  sln0SS 
Zqq  =  -nig  COS  0-q 


$ESS=  USS3in4'sSc°sdSS  +  V 

tan0ss  =  wss^uss 


ss 


cos4,SS+WSSsin4,SSsiri9SS- 


(u 


ss 


w  ^ ^ 

oS  ' 


Since  bank  angle  equals  zero  and  the  sum  of  flight 
path  angle  and  angle  of  attack  equal  the  pitch  angle,  then 
in  steady  state  pitch  angle  equals  angle  of  attack. 

0  =  <X 

SS  ss 

These  equations  are  used  to  define  the  initial  cruise  condi¬ 
tion  of  the  missile  which  is  low  level  straight  and  level 
cruise  flight  at  0.75  Mach  number.  The  precise  initial 
flight  condition  and  initial  target  position  is  described  in 
Appendix  B. 

5.  Other  Useful  Relations 

Components  of  load  factor  at  the  missile  center  of 
mass  in  the  y  and  -z  directions  are  given  below: 


n 


z 


(10a) 


V 


(S  +  R) 


y  <ng  g 

6.  Definition  of  Controls  and  Control  Limits 

The  control  configuration  considered  in  the  thesis 


(10b) 


is  a  conventional  wing  and  tail  with  conventional  rudder 
(0,  aileron  (§)  and  elevator  (h)  control  surfaces.  In 
addition,  to  cover  a  missile  configuration  which  uses 


differential  horizontal  stabilizer  (stabilator)  for  roll 
control  as  opposed  to  aileron,  the  quantities  left  stabilator 
deflection  (n^)  and  right  stabilator  deflection  (n^)  are 
defined  and  calculated. 

Limits  are  imposed  on  rudder,  elevator,  and  aileron 
travel.  Limits  are  not  imposed  on  differential  stabilator 
travel  but  stabilator  travel  is  presented  for  all  simulations 
to  show  the  stabilator  travel  required  if  indeed  that  were 
the  control  configuration.  The  aileron  power  coefficient 
data  presented  in  Appendix  A  reflects  the  limited  control 
authority  in  roll  typical  of  a  missile  configured  with 
differential  stabilator  for  roll  control  as  opposed  to 
aileron. 

The  following  control  definitions  and  limits  are 

applied : 

a)  Thrust  (T)  is  a  constant  and  is  oriented  along 
the  x  axis. 

T  =  T  =  constant 
x 

b)  Rudder  deflection  (s)  is  defined  as  positive 
trailing  edge  left  from  the  rear  and  is  limited  to 

-15°  <  C  £  15° 

c)  Elevator  deflection  (n)  is  defined  as  positive 
trailing  edge  down  (produces  a  nose  down  moment)  and  is 

-15°  <  n  <  15° 


limited  to 


d)  Aileron  deflection  (£.)  is  defined  as  positive 
with  trailing  edge  of  the  right  aileron  up  (produces  a 
positive  (clockwise  from  the  rear)  roll)  and  is  limited  to 

-15  <■  5  <  15° 

e)  Left  horizontal  stabilizer  deflection  ( )  is 
positive  trailing  edge  down. 

nL  =  n  +  C 

f)  Right  horizontal  stabilizer  deflection  (n_)  is 
positive  trailing  edge  down. 

nR  =  n  -  S 

B.  AERODYNAMIC  COEFFICIENT  BUILD  UP 


Aerodynamic  forces  and  moments  are  built  up  from  coeffi¬ 
cient  data  in  standard  fashion  using  the  following  relations: 
D=C^qS  L  =  CT  q  S 

M.  -  CM  qSc 


Y  =■  Cy  5  S 


‘A  ''M  = 


NA  =  Cn  q  Sb 


1  •  Definition  of  Coefficients 

Each  coefficient  is  partitioned  into  two  parts;  a 
static  term  and  a  dynamic  term  as  follows: 

u  UST  UDYN 

CT  =  CT  +  CT 

L  bST  LDYN 

Dy  =  Cy  +  Gy 

1  oT  1DYN 

c  =  cm  +  cm 

m  mST  mDYN 

C.  =  C,  +  C-, 

1  ■lST  •lDYN 
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2.  Coefficient  3uild  Up 

The  static  terms  are  constructed  from  coefficient 
data  stored  in  tables  (Appendix  A)  as  functions  of  two  or 
three  variables  (either  «,  3,  M,  n»  C*C)  as  follows: 


Jbasic 


=  f  (  *,  8.M) 


ACt  (d) 

lst 

=  f 

( 3*  n) 

Dbasic 

f  ( 

<r,  g,M) 

ACn  (n) 
UST 

=  f 

( <*»  0,  n) 

ACr,  (C) 
UST 

=  f 

acd  (c) 
UST 

=  f 

(**8*5) 

CY  = 

basic 

f  ( 

oc,  3) 

acy  (^) 

X3T 

=  f 

( «»  B*  5) 

ACV  (c) 
*ST 

=  f 

( 3*  5) 

c 

mbasic 

f  ( 

«» 3t M) 

aC  (n)  =  f  («.s»n) 


ST 


n 


basic 


=  f  (  <*>  8) 


AC  (5)  =  f  («,3»q) 
nST 


AC  (5)  =  f  («*8»5) 


=  f  («,6) 


'basic 


C,  (?)  =  f  («.s,e) 

XST 

C,  (?)  =  f  («,8,?) 

XST 

The  static  coefficients  are  formed  as: 


JST 


:T  +  act  (n) 

L basic  L 


=  CT 


ST 


Cy  =  Cy 

XST  1  basic 


m 


=  C 


ST 


m 


basic 


+  AC-  (n)  +  AC 

(  O  +  AC 

ST  UST 

+  ACV  (?)+ACv 

(€) 

iST  iST 

+  ac  (n) 

ni  Q  m 
b  i 

D 

JST 


C  n  =  C  n  +  AC  n  (  C)  +  ACn  (  ?) 

ST  basic  ST  ST 

C,  =  C,  +  AC,  (S)  +  AC,  (?) 

XST  xbasic  X3T  XST 


The  dynamic  terms  are  constructed  from  dynamic 
coefficient  data  stored  in  tables  (Appendix  A)  as  a  function 
of  angle  of  attack  and  Mach  number  f  (®,M).  They  are 
ct  .*  ct  »  cn*»  cn  *  cv.»  Cy  ,  Cy  ,  Cm  C,  ,  cnv  cn  »  C  r,  ' 


T  •» 

L  Lq 


'D  WD  ’  Y 


8  xr  xp  moc  ^q 


C-^  ,  and  • 


The  dynamic  coefficients  are  formed  as: 

ct  =  —  (CL.«  +  C  q) 

LDYN  2Vt  L «  q 


DYN  2Vr 


( CD i"  +  CD  q) 


'DYN  2V 


(  Cy  r  +  Cy  P  +  Cy  •  6) 
*  Ya 
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n 


DYN 

=  -2-  (G  «  + 
2Vt 

C 

mq 

q) 

=  —  ( G  3  + 

G 

r  +  C 

DYN 

2Vt  "6 

nr 

=  —  (C,  r 

+  C, 

P) 

DYN 

2Vt  xr 

X 

P 

n 


The  total  coefficient  build  up  is  summarized  below. 


Ct  =  1  “UL 

L  L basic  LST 


+  ACt  (  n)+  -  (CT  . i+CT  q) 

L“  Lq 


2Vr 


+ACV  (c)+ACy  (£)  +  —  (C  r+Cy  p+Cy.S) 


-'v  —  w  V  '  —  V 

1  Y basic  ST 


'ST 

c" 


2Vt  "  xr  3 


C  =  C  +AC  (n)+  — -  (C  «+C  q) 

31  m  basic  ST  2Vt  m«  mq 

C  =  C  +AC  (C)+AC  (5)+  2?“(C  §  +  C  r+C  p) 

n  ncasic  nST  nST  2VT  n§  nr  np 


C,  =  C,  +AC,  (c)+aC,  (g)+  — 

1  1 basic  XST  '  X3T  2V, 


■(C,  r+C.,  p) 
r  ~d 


C.  MISSILE  LINEARIZED  EQUATIONS  OF  MOTION 

The  linearized  aquations  of  motion  used  to  design  the 
missile  autopilot  are  presented  in  this  section.  The  equations 
are  linearized  about  the  cruise  flight  condition. 

1 .  Additional  Assumptions  for  Linearized  Equations 

The  additional  assumptions  required  for  linearized 
equations  are  as  follows: 

a)  Derivatives  are  given  in  stability  axes  which 
are  fixed  for  the  missile  low  level  cruise  condition. 

b)  Standard  small  perturbation  assumptions  apply. 

1)  Perturbed  angles  are  small. 


JT. 
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2)  Products  of  perturbed  variables  are  neglected. 

3)  There  is  no  coupling  between  longitudinal 
and  lateral  directional  motions. 

c)  Thrust  derivatives  are  neglected. 

d)  The  steady  state  condition  is  chosen  as  low 
level  straight  and  level  cruise. 

2.  Linearized  Equations  Summarized 

The  linear  equations  of  motion  in  state  variable 
form  are  shown  in  Tables  2-1  and  2-2.  The  following 
coefficients  are  defined  as  used  in  Table  2-2. 


A 


^ZZ^XZ 


I  I  -I 

ZZ  XX  xz 


c  = 


^xx^xz 


^zz^xx~^xz 


3.  Definition  of  Stability  Derivatives 

The  dimensional  aerodynamic  stability  derivatives 
used  in  the  equations  of  motion  are  defined  by  the  following 
relations. 


a.  Longitudinal  Dimensional  Stability  Derivatives 


xu 

'qSSS(CDl7+2CDss)/mlJSS 

1  /sec 

3  q*SSS(CLss"CDoc)/in 

ft/ sec 

X 

n 

3  -qssSCD  /a 
n 

ft/sec 

ZU 

3  qSSS(CLu+2CLss/mUSS 

1  /sec 

33 


Longitudinal  Equations  (Perturbed  Variables)  (11) 


Z<*  =  "qSSS(CL  +CDss^/m 

ft/sec2 

Z<S  =  ■qSSS5>CL.//2mUSS 

ac 

ft/sec 

Zq  =  "qSSS°CL  ^2mUSS 

q 

ft/ sec 

Zn  =  “qssscLn//m 

ft/sec 

MU  =  qSSSc(CMu  +  2CMss)/IyylISS 

1/ft  sec 

M«  "  *SSSBCM«/Iyy 

1 /sec2 

»i  -  ?SSS52<V2IyyUSS 

1  /sec 

Mq  ■  5sSSj2c»q/2IyyII5S 

1  /sec 

"n  ■  ’SS^'^/^y 

1  /sec 

b.  Lateral  Directional  Dimensional 

Derivatives 

Stability 

YS  =  5SSSCyB/m 

ft/ sec2 

Yp  ”  ^SSSbCy  ^2aUSS 

P 

ft/sec 

Yr  =  qSSSbGy  //2mUSS 
r 

ft/ sec 

YS  =  W°y/S 

ft/ sec2 

YC  =  «SSSCy/" 

ft/ sec2 

L0  =  ^SSSb0l8/Ixx 

1 / sec2 

lp  1  *SSSb2ci  /2IxxJSS 

P 

1  /sec 

tr  -  4SSSb2ci  /2IxxDSS 

1  /sec 
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1  /sec 


h  =  qSSSbClc/rxx 
h  *  ?SSSb0l/Ixx  1/sse2 

h  *  ^SSSh\nzz  1/se=2 

%  *  7SSSb\  /2IzxDSS  1/3eo 

P 

“r  *  5ssSb2Cn  /2IZVSS  1/«c 

r 

N€  =  7ss3bCns/I2z  1 /sec2 

\  =  «SSSMn 

s 


4.  Trim  Equations 

For  straight  and  level  flight  at  constant  velocity 
the  following  variables  are  all  zero:  U  =  V=W=p=Q  =  R=9  =  '}<=$ 


=ZS  =  0  QSS=PSS  RSS  =  <I,SS=LA3S=MA33=NAS3  yss  yss  =  bss  =  vss: 


=  QC  =  0 


ss  ^ss  Sss 

The  trim  equations  reduce  to: 


Tsin“sS  -  «g-(0L  +CL  -sstCL  nss)7ssS 


TC0S“SS  "  (CD  +CD _rSS+CDnnSS,^SSS 
o  «  n 


°  =  (V+0m  _*SS+C.  ’ ^SS>^SSS~ 


r 


III.  DEVELOPMENT  AND  SIMULATION  OF  THE  MISSILE  AUTOPILOT 

The  missile  is  assumed  to  be  equipped  with  an  autopilot 
capable  of  providing  closed  loop  control  which  consists  of 
the  following: 

a)  normal  acceleration  (n  ) 

z 

b)  lateral  acceleration  (n  ) 

y 

c)  bank  angle  ( <t>) . 

In  addition,  an  outer  closed  loop  control,  which  serves 
as  the  outer  autopilot,  is  as  follows: 

a)  altitude  (H) 

b)  vertical  flight  path  angle  (y) 

These  two  outer  loops  are  employed  as  required  during  certain 
phases  of  the  attack  mission. 

It  is  assumed  that  the  missile  possesses  accelerometer 
sensors  in  the  y  and  z  body  axes  located  at  the  missile 
center  of  gravity,  rate  gyros  and  rate  integrating  gyros 
about  the  x,  y,  and  z  axes,  and  a  radar  altimeter.  Sensor 
noise  is  neglected  in  this  simulation. 

A.  AUTOPILOT  INNER  LOOP  REQUIREMENTS  AND  DESIGNS 
1 .  Description  of  the  Autopilot  Inner  Loops 

The  autopilot  is  assumed  to  employ  three  inner  loops 
as  depicted  in  Table  3-1. 
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TABLE  3-1 

Autopilot  Inner  Loop  Description 


LOOP 

COMMANDED 

CONTROL 

CONSTRAINTS 

VARIABLE 

1  Normal 

n  (normal 

n 

-15°  * 

n  *  15° 

Acceleration 

z 

Command 

acceleration ) 

El evator 

-2  S 

N,  -  4 

4J 

2  Bank 

C 

1 

— k 

VJ1 

o 

IA 

5  -  15° 

Angle 

Command 

(Bank  Angle) 

Aileron 

4>  S 

60°* 

p  ^  50,1 00,200°/sec 


3  Turn 

n  (Lateral 

4 

Coordinator 

y 

Acceleration) 

Rudder 

-15°  -  ?  ^  15° 

N  =  0 

y 


*  Applies  in  certain  mission  phases  only 
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2.  Normal  Acceleration  Command  System 

The  purpose  of  normal  acceleration  command  (NAC) 

system  is  to  provide  a  vertical  load  factor  (n2)  response  to 

a  commanded  load  factor  (n  ).  An  accelerometer  is  used  to 

zc 

provide  the  primary  feedback  and  a  rate  gyro  is  used  to 
provide  inner  loop  pitch  rate  feedback  for  improved  damping. 

a.  Block  Diagram 

The  block  diagram  of  the  normal  acceleration 
command  autopilot  is  shown  in  Figure  3-1 .  Limiter  a  allows 
the  commanded  acceleration  to  a  range  from  -2.0  to  +4.0  g’s. 
Limiter  b  allows  the  elevator  control  to  vary  from  -1 5 
degrees  to  +15  degrees. 

b.  Assumptions 

1 )  The  rate  gyro  and  accelerometer  dynamic  lags 
are  negligible. 

2)  The  accelerometer  is  mounted  at  the  center 
of  gravity.  Therefore  the  moment  arm  is  zero  and  the 
feedback  load  factor  is  totally  nz*  (C=0) 

3)  The  fix  servo  can  be  represented  by  a  first 

order  lag. 

4)  The  missile  dynamics  can  be  represented  by 
a  short  period  approximation. 

5)  Commanded  normal  acceleration  is  limited  from 
-2.0  to  +4.0  g's. 

c.  Design 

The  following  transfer  function  was  developed 
from  equation  (11): 
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VyZnMi)»  *  (M.Zn-Z.Mn) 

VT{  s2-  ( M  +Z  +  M-)  s  +  (Z.M  -Mj  i 
i  q  «  oc  q  J 

VT  VT 

From  this  transfer  function  and  equation  (10a) 


n 

z 


(Q-« ) 


the  following  9  to  nz  transfer  function  was  obtained. 


V£ 

Q(s) 


g  ( VmM  +Z  M. )s+(M  Z  -Z  M  ) 

if)  n  «  crqocfj 


Utilizing  the  characteristics  listed  in  Appendix  A 
for  this  cruise  missile  at  0.75  Mach  number,  the  above  transfer 
functions  become 


Q ( s )  =  -26.99(3+1.009)  and  N2(S)  =  -0.08(s2+0.1985s-323.9) 
n(s)  s2+1 .289S+27.1 3  Q(s)  (s+1.009) 

Evaluating  the  characteristic  polynomial  of 
the  Q(s)/p(s)  transfer  function  yields  the  natural  frequency 
of  0.829  Hz  and  a  damping  ratio  of  0.1237.  The  resulting 
block  diagram  for  the  normal  acceleration  command  inner 
loop  autopilot  for  analysis  purposes  is  shown  in  Figure  3-2. 

( 1 )  Root  Locus  Evaluation. 

(a)  Pitch  Rate  Loop.  The  transfer  function 
for  the  pitch  rate  loop  is: 
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Q(s)  1080  (s-H) _ 

E(s)  (s3+41 .29s2+(78.7+1080KR1 ) s+ ( 1  084+ 1  080KR1  ) 

Appendix  C  contains  the  root  locus  plot  of 
the  pitch  rate  loop  for  ranging  KR1 .  A  KR1  of  0.28  is 
chosen  which  yields  the  following  characteristics: 
s  =  -5.600  ±3.844  j 
<t>  =  55.  53° 

5  =  0.82 

u>  =  1  .08  H 
n  2 

(b)  Normal  Acceleration  Command  Loop. 

Using  a  KR1  of  0.28  in  the  pitch  rate  loop*  and  locating  the 
acceleration  at  the  missile  center  of  mass  (C=0),  the  normal 
acceleration  command  loop  becomes: 


-86.4KA1 (s2+0.2s-329) 

f  N 

_ _3 _ V 

s(s3+41  . 29s  2+  381  ,1s-1  386.4) 

L _ 

N  (s)  -86.4KA1 (s2+0.2s-329) 

N  (s)  s(s3+41 .29s+381 .1 s-1 386.4)-86.4KA1 (s2+0.2s-329) 

zc 

Appendix  C  contains  the  root  locus  plot  of 
the  normal  acceleration  command  system  for  varying  KA1 .  A 
KA1  of  0.05  is  chosen  which  yields  the  following  charac¬ 
teristics: 

s  *  -1.981  ±  3.317  j 

t>  =  30.85° 

5  =  0.513 

=  0.61  5H 

n  z 
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With  a  KR1  of  0.28  and  a  KA1  of  0.05  the 


system  transfer  function  relating  n  to  n  is: 

zc  z 

N  ( s)  -4.32(s2-0.2s-329) 

z _  _ 

N  (s)  s4+41.29s3+376.78  s2-1  387.26s-U21 .28 

z 

c 

The  normal  acceleration  command  autopilot 
has  zero  steady  state  error  for  a  step  input  as  it  is  a 
type  1  system.  The  system  response  to  a  1 g  commanded  load 
factor  is  shown  in  Figure  3-3. 

3.  Bank  Angle  Command  System 

The  purpose  of  the  bank  angle  command  (BAC)  auto¬ 
pilot  is  to  command  missile  bank  angle.  Limits  are  applied 
to  commanded  roll  rate  and  in  certain  mission  phases 
commanded  bank  angle.  A  rate  integrating  gyro  feeds  back 
bank  angle  and  a  roll  rate  gyro  provides  feedback  as  a 
roll  rate  damper. 

a.  Block  Diagram 

The  bank  angle  command  autopilot  is  presented 
in  block  diagram  in  Figure  3-4.  Limiter  a  allows  the 
commanded  bank  angle  to  vary  from  0  to  60  degrees  and  is 
active  in  certain  flight  phases.  Limiter  b  limits  the  roll 
rate  to  either  50,  100,  or  200  degrees/second.  Limiter 
c  limits  aileron  deflections  to  range  from  -15  to  +15 
degrees . 


b.  Assumptions 

1  )  The  rate  gyro  and  displacement  dynamic  lags 


are  negligible. 
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2)  The  fin  servo  can  be  represented  by  a  first. 


order  lag. 


3)  The  missile  dynamics  can  be  represented  by 


a  roll  mode  approximation. 

4)  There  are  limits  on  commanded  bank  angle  and 
commanded  roll  rate  as  noted  above, 
c.  Design 

Equation  (12)  can  be  written  to  find  the 
lateral  transfer  functions  ^  |  ^  j-,  \  *  and  ~~j*  * 

where  5  =  £  or  5.  By  allowing  a  roll  approximation  the 

transfer  function  relating  aileron  and  bank  angle  becomes: 

<J)  ( s)  L  ^ 

5(s)  s  (  s-L  ) 

F 

The  root  at  the  origin  indicates  the  missile  does  not  have 

bank  angle  stability  without  autopilot  stabilization. 

Evaluation  of  the  transfer  function  coefficients  results 

in  the  transfer  functions  ~ ^ ^  and  ^  s  ^  = 

E, ( s)  s(s+6.2)  C 


-116 
s+6. 2 


The  maximum  missile  roll  rate  can  be  approximated 


:  ii6)(i 


=  280.6  degrees/second. 


The  block  diagram  of  the  bank  angle  command 

system  is: 
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(1  )  Root  Locus  Evaluation. 

(a)  Roll  Rate  Loop.  The  transfer  function 
for  the  roll  rate  loop  is: 


P(s)  _  1640  KR2 _ 

E(s)  s2+46.2s+248+4640KR2 


Appendix  C  contains  the  root  locus  plot 
of  the  roll  rate  loop  for  varying  KR2 .  A  gain  KR2  of  0.10 
yields  the  following  characteristics. 

s  =  -23.10  ±  13.36  s 
<J>  =  59.96° 

£  =  0.8656 

ui  =  4.247H7 
n  L 

(b)  Bank  Angle  Command  Loop.  With  a  02 
of  0.10  the  bank  angle  command  loop  becomes: 


__t r\ _ 

464KD 

4> 

_ s. 

s 

4 

\ 

s(s2+46.2s+712) 

- -7 

The  transfer  function  is: 

<fr(s)  _  464KD _ 

6  (s)  s3+46.2s5+712s+464KD 

c 

The  root  locus  plot  is  shown  in  Appendix  C 
for  varying  KD.  Choosing  a  KD  gain  of  10.8  yields  the 
following  characteristics: 
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■9.98  ±  9 . 5 8 j 


$  =  45.056° 
C  =  0.718 
4.  =  2.19  H„ 


Using  gains  of  0.10  for  the  inner  loop  and 
10.8  for  the  outer  loop  of  the  bank  angle  command  system, 
the  system  transfer  function  G{s)  becomes 


*cU) 


5011  .2 


s'5  +  46.2s‘i+712s+5011  .2 


and  the  system  is  a  type  1  system  for  a  step  input.  The 
transient  reponse  is  demonstrated  in  Figures  3-5  and  3-6 
for  a  commanded  bank  angle  of  60  degrees. 

4.  Turn  Coordinator 

The  turn  coordinator  (TC)  is  a  lateral  autopilot 

which  provides  a  body  directional  load  factor,  n  ,  to  a 

commanded  load  factor,  n  .  Like  the  normal  acceleration 

^c 

command  autopilot  it  employs  an  inner  loop  yaw  rate  damper. 
The  outer  loop  uses  unity  accelerometer  feedback  and  has 
the  capability  of  moment  arm  feedback  if  the  accelerometer 
is  located  at  a  different  position  than  the  missile  center 
of  mass. 

a.  Block  Diagram 

The  block  diagram  representing  the  turn  coordina¬ 
tor  autopilot  is  shown  in  Figure  3-7.  Limiter  a  limits 
the  rudder  deflection  (^)  from  -15  to  +15  degrees. 


r  f s) 

C  <s>  VTs2-YetVTNr)s+(YsNrtV,tNs.YrN8) 

Using  the  equation  (10b)  in  Chapter  I  the 
transfer  function  for  yaw  rate  to  n  is 

Ny(a)  _  Y,s2^(YrNr-YrNr)s+(I^NB-Y3Nc) 

R(S)  g  VTV+(YrN3-YS\) 

Substituting  coefficient  data  obtained  from  Appendix  A 

results  in  the  transfer  functions  as: 

R ( s )  =  -1 4. 5(s+0.11 ) 

Z ( s )  s  +0 . 245s+1 4-5 


N  (s)  -0.1 4(s2-0.01 s-20.9) 

R  (s)  s+0.11 

The  resulting  turn  coordinator  block  diagram 
for  analysis  purposes  is  given  in  Figure  3-8. 

( 1  )  Root  Locus  Evaluation. 

(a)  Yaw  Rate  Loop.  The  loop  transfer  func¬ 


tion  is 


R(s)  580(s+0.11  ) _ 

E ( s )  s3+40.25s*+(24.3+580KR3)s+(580+63.8KR3) 

Appendix  C  shows  the  root  locus  plot  of 
the  yaw  rate  loop  for  varying  KR3.  Selecting  a  gain  KR3 
of  0.40  results  in  the  following  characteristics: 

s  =  -3.581  ±  2. 339j 
<j>  =  56.85° 

,  =  0.8372 
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KA2 


W  *  0.6807H 
n  z 

<j>  =  56,85  degrees 

(b)  Turn  Coordinator  Loop.  With  a  KR3  of 
0.40  the  turn  coordinator  loop  becomes: 


-81 .2KA2(s2-0.01 s-20.9) 

*  - 

r  | 

s(s^+40.25s^+256.3s+605.52) 

— 3 - 7 

N  (s) 

I— 

N  (s) 
^0 


The  transfer  function  is: 

-81 ,2KA2(s2-0.01 s-20.9) 

s4+40.25s3+(244.7-81 . 2KA2 ) s2+ ( 604.2+81 ,2KA2)s  + 

1697KA2 


The  root  locus  plot  for  varying  KA2  is  shown 
in  Appendix  C.  Selecting  a  gain  KA2  of  0.25  yields  character 
istics  as  follows: 

3  =  -2.206±  2.807J 
4>  =  38.16° 

5  =  0 . 6l 8 

ujn  =  0.5682Hz 

The  turn  coordinator  system  transfer  function 

G(s)  becomes: 

Ny(s)  -32. 48(s2-0. 01 s-20.9) 

N  (s)  s4+40.25s3+21  2.2s2+604.  5s+678.8 

yc 

The  system  is  type  1  and  therefore  has  zero  steady  state 
error  to  a  step  input.  Figure  2-5a  shows  the  response 
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of  the  turn  coordinator  for  a  step  aileron  input.  Note  that 


an  n  of  zero  results  in  a  coordinated  turn. 

^c 

B.  AUTOPILOT  OUTER  LOOP  REQUIREMENTS  AND  DESIGNS 

1 .  Description 

The  outer  loop  auto  pilot  has  two  functional  loops 
consisting  of  an  altitude  hold  and  a  vertical  flight  path 
angle  hold. 

2.  Vertical  Flight  Path  Angle  System 

The  purpose  of  the  vertical  flight  path  angle  hold 
loop  or  gamma  command  loop  is  to  provide  missile  reponse  to 
a  command  of  8.6  degrees  during  a  specified  phase  in  the 
terminal  profile. 


a.  Block  Diagram 


b.  Assumption 

1 )  Bank  angle  is  zero. 

c.  Design 

Starting  with  equation  (10a)  n  =  T(Q-g)  =  _^T  y 

g  g 

the  transfer  function  =  V^Ts)  can  be  deriveci*  Using 

z  T 

the  inner  loop  design  previously  for  load  factor  the  vertical 
flight  path  angle  block  diagram  becomes: 


1_ 

s 


v  .n 

R 

-9.504(s2+0.2s-329) 

c  1 

1 

s4+41 .29s^+371  .6s2+1384s+3127 

Y  (  S  ) 

The  transfer  function  “yTs)  ^s: 

-9.504K  (s2+0.2s-329) _ _ 

s5+41.29s4+37.6s-i+(  1384-9. 504K  ) sS+  (31  27-1  .  9K  ) s+31 27K 


(1)  Root-Locus  Evaluation.  Appendix  C  shows  the 
root  locus  plot  of  the  vertical  flight  path  angle  loop  for 
varying  gain  K  .  A  gain  Ky  of  1.0  yields  the  following 
characteristics : 

s  =  1 .2725  ± 2. 5407 j 

(D  =  26.6° 

;  =  0.4478 

un  =  °‘45HZ 

The  system  transfer  function  l(s)  with  a 


gain  of  1  . 0 

yL5),  = 
Yc  (a) 


becomes: 

-9. 504(s2+0.2s-329) _ _ _ 

s5+41  .29s4+371 . 6s ^+1  374.  5s 2+ 31  25 . 1  s+31  27 


Figure  3-9  depicts  the  system  response  to 
a  commanded  vertical  flight  path  angle  of  3.6  degrees. 

3.  Altitude  Hold  System 

The  purpose  of  the  altitude  hold  autopilot  is  to 
hold  altitude  at  50  feet  during  specified  flight  phases, 
a.  Block  Diagram 


b.  Assumptions 

1)  Vertical  rate  of  climb  (or  desent) 
subtends  a  small  angle  (  y is  less  than  20°),  therefore 


H  =  VTsiny  =  V?Y . 

2)  Bank  angle  equals  zero, 
c.  Design 

Let  altitude  command  and  altitude  be  designated 
and  H  respectfully.  The  design  is  based  upon  the 
following  two  relations: 


H  =  Kh(Hc-H) 

'  *  in  (H„-H) 


From  the  vertical  flight  path  angle  development, 

Y(s)  _  g _ 

Nz(s)  VTs 

sH(s)  =  VTy(s)  (Divide  both  sides  by  an(s)) 


sH(s ) 

_  -  V-, 

Y  (_S  )  =  VT  _g _ 

Va> 

T 

N^Ts)  T  VTs 

a(s) 

II 

Y  (s) 

s 

The  altitude  hold  block  diagram  becomes: 


The  transfer  function  H(s)/Kc(s)  is: 


HU)  = 
HCU) 


-9. 504KH( s  tQ.2s-329 ) _ 

s6+41 .29a^+371 .6s4+1374.5s3+(3125.1-9.504KH 
(3127-1 .9KH)s+3127KH 


(1 )  Root  Locus  Evaluation.  Appendix  C  contains 
the  root  locus  plot  of  the  altitude  hold  loop  for  a  varying 
gain  KH.  Again  KH  of  0.3  results  in  characteristics: 
s  =  -1.4045  ±  2.424j 
4>  =  30.08° 
e  =  0.501 
wn  =  0.445  Hz 

The  altitude  hold  loop  transfer  function 
G(s)  with  KH  of  0.3  is: 

H(a)  =  2.8512(s2+0.2s-329) _ _ _ 

Hc(s)  s>  41 .29s  ^+3  71  . 6s4+1  374- 5s3+31  22. 2  5s2+31 26. 4s+938 . 1 


and  Figure  3-10  shows  the  response  to  a  y  c  of  50  feet. 
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IV.  DEVELOPMENT  AND  SIMULATION  OF  MISSILE  GUIDANCE  SYSTEM 

The  purpose  of  this  study  is  to  evaluate  the  effect  of 
limited  roll  rate  and  various  guidance  laws  on  the  terminal 
performance  of  a  bank-to-turn  cruise  missile  using  propor¬ 
tional  navigation  in  the  terminal  (attack)  phase  against  a 
medium  to  large  combatant  ship.  A  baseline  guidance  and 
control  scheme  is  designed  and  studied  which  uses  a  pop 
out  maneuver,  proportional  navigation,  and  bank-to-turn 
maneuvering  in  the  attack  phase.  The  trajectory  parameters 
and  navigation  constants  are  optimized  for  best  performance 
against  a  moving  ship  utilizing  blinking  countermeasures. 

The  missile  is  assumed  to  have  limited  roll  performance 
due  to  the  use  of  differential  stabilizer  for  roll  control 
vice  ailerons.  The  terminal  performance  of  this  baseline 
system  is  evaluated  for  varying  roll  rate  limits,  electronic 
countermeasure  blinking  rates  and  burn  through  ranges. 

A  second  control  scheme  is  also  designed  and  evaluated 
for  the  attack  phase.  This  mode  uses  a  sea  skimming  scheme 
in  which  a  50  foot  altitude  is  maintained  and  proportional 
navigation  is  used  in  azimuth  during  ingress  until  a  specified 
range  is  reached  at  which  time  proportional  navigation  is 
used  in  both  azimuth  and  elevation  to  attack  the  target. 
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A.  MISSION  DESCRIPTION 

1 .  Baseline  Guidance  Scheme  (Pop  Out  Maneuver) 

The  mission  is  divided  into  four  phases: 

a)  ingress 

b)  offset 

c)  popup 

d)  attack 

The  mission  commences  with  target  ship  located  25,000  feet 
north  of  the  missile  and  moving  east  at  21  knots.  The 
missile  initial  heading  is  north. 

'ihe  ingress  phase  commences  with  the  missile  flying 
at  an  altitude  of  50  feet  on  altitude  hold  and  utilizing 
proportional  navigation  in  azimuth  to  home  in  on  the 
target.  At  a  given  range  the  offset  phase  is  initiated  with 
a  turn  to  the  right  of  approximately  12  degrees  followed  by 
a  commanded  bank  angle  of  zero  degrees.  No  proportional 
navigation  is  used  in  this  phase.  At  another  given  range 
a  pop  up  maneuver  is  initiated  which  generates  a  low  angle 
climb  in  preparation  for  the  attack  phase.  During  the  pop  up 
maneuver,  proportional  navigation  is  used  to  begin  turning 
the  missile  toward  the  target.  Finally  the  missile  enters 
the  attack  phase  at  a  specified  range  where  proportional 
navigation  is  used  in  both  vertical  and  azimuth  planes. 

In  this  phase  the  bank  angle  is  unlimited  but  roll  rates  are 
limited  and  used  as  one  of  the  parameters  in  studying  miss 
distances.  The  second  parameter  under  study  is  the  electronic 
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counter  measure  blinking  rate  frequency  which  is  initiated 
at  the  start  at  the  start  of  the  attack  phase  and  allowed  to 
remain  activated  until  a  specified  burn  through  range  which 
is  the  third  miss  distance  study  parameter.  Figure  4-1  is 
a  geometric  depiction  of  the  missile  flight  path. 

2.  Alternate  Guidance  Scheme  (Sea  Skimmer) 

The  mission  for  this  guidance  scheme  is  divided  into 
two  phases: 

a.  ingress 

b.  attack 

This  mission  begins  with  the  same  initial  scenario  as  the 
baseline  scheme,  however,  the  pop  out  maneuvers  (turn  and 
pop  up  phases)  is  eliminated. 

The  missile  flies  at  50  foot  altitude  during  the 
ingress  phase  and  utilizes  proportional  navigation  in 
azimuth  until  a  range  of  3000  feet  at  which  point  the 
missile  enters  the  attack  phase  and  dives  in  on  the  target 
using  proportional  navigation  in  both  azimuth  and  elevation. 
This  low  altitude  approach  essentially  resembles  a  sea 
skimming  effect.  Figure  4-2  shows  the  geometric  flight 
profile  of  this  scheme. 

3.  Electronic  Counter  Measures  (ECM) 

It  is  assumed  that  the  target  has  the  ECM  capability 
of  simultaneously  shifting  the  aim  point  of  the  missile  75 
feet  aft  and  10  feet  vertically.  The  simulation  is  tested 
for  ECM  blinking  frequencies  of  0.5,  1.0,  and  2.0  cycles 


k 
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per  second.  Blinking  is  commenced  during  the  attack  phase 
and  continues  until  jammer  burn  through  occurs  when  the 
target  signal  power  can  be  seen  over  jammer  power.  Although 
3800  feet  (as  calculated  using  Hosington,  {Ref.  3})  is 
considered  a  typical  burn  through  range  the  missile  is 
tested  for  burn  through  ranges  from  3800  feet  down  to 
400  feet. 

4.  Glint 

This  study  models  glint  with  a  gaussian  function 
that  gives  the  aim  point  a  random  shift  effect  with  a 
variance  of  25  feet  in  the  Y  axis  and  4  feet  in  both  the 
X  and  Z  axes.  Since  glint  is  imposed  as  a  fluctuation  in 
target  coordinates  and  target  coordinates  are  used  to 
calculate  the  target  line  of  sight  angles  then  glint  is 
modeled  as  inversely  proportional  to  range  as  required. 

For  ease  in  simulation  glint  is  imposed  only  during  the 
attack  phase.  Results  for  runs  with  glint  included  are 
presented  in  Chapter  IV. 


3.  SEEKER  EQUATIONS  AND  SIMULATION 

A  flow  diagram  of  the  seeker  simulation  is  presented 
in  Figure  4-3. 

1 .  Line  of  Sight  Rates 

It  is  assumed  that  the  missile  seeker  has  the  ability 

to  sense  target  range  and  target  line  of  sight  rates  in 

azimuth  and  elevation  ( and  'the  body  coordinate 

B  B 

system  out  to  a  range  of  2 5,000  feet. 
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These  line  of  sight  rates  are  generated  in  simulation 
by  a  two  step  procedure.  First,  the  line  of  sight  rates 
in  earth  referenced  azimuth  and  elevation  (  a  ^  and  a 
are  constructed  by  using  the  following  relations: 


^  T  =  tan 


-1  XT  -  YE 
XT  “  XE 


9_  =  tan"1  — - - p - 5— 

1  ((xT-xEr+(YT-YEr) 

O  _  ur  _  u» 

AZ  T 

•  •  • 

"AZ  '  V* 

°EL=  V9 

S  =  S  .0 

EL  T 


Figures  4-4  and  4-5  depict  the  geometry  represented  by 
the  above  equations. 

Second,  the  body  line  of  sight  rates  are  calculated  from 

the  following  transformation. 

•  •  • 

oAZ  =  -  aEEsin<J>+  aAZ  cos  9cos  <i> 

5 

•  •  • 

aEL  =  a^2cos9s^n<^ 

B 

.  •  •  . 

Once  the  body  rates  have  been  simulated  ( o ^  ,  ct,^  ) 

B  B 

they  are  transformed  to  earth  horizontal  and  vertical  rates 
through  the  transformation  T defined  by: 

'■el  ■  %LBeos*-5AZBsin* 


b  .  b 

ft  7.  ITT 


sin<to-o,„  cos<t> 


a 


AS 


Figure  4-4  Azimuth  Line  of  Sight  Description 
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This  transformation  includes  the  assumption  that©  is  small 
throughout  the  mission.  The  earth  horizontal  and  vertical 
line  of  sight  rates  (a^  and  a^)  are  then  used  in  the 
appropriate  guidance  laws  for  proportional  navigation. 

2.  Missile  Dynamic  Filters 

In  order  to  filter  out  missile  dynamics  from  the 
seeker  line  of  sight  rates,  a  low  pass  filter,  10/(S+10), 
is  used  at  the  input  to  the  seeker  on  0  and  $ . 


C.  BASELINE  GUIDANCE  LAW  DESIGN  (POPOUT  MANEUVER) 

A  flow  diagram  of  the  baseline  guidance  scheme  is 

shown  in  Figure  4-6.  'Switch  positions  are  a  function  of 

\ 

mission  phase.  Each  phase  of  the  mission  is  described 
below  beginning  with  the  attack  phase. 

1.  Attack  Phase 


Proportional  navigation  is  used  in  azimuth  and 
elevation  to  construct  guidance  commands  to  the  missile 
autopilot  commencing  at  a  range  of  9100  feet  in  the  manner 
described  below. 

a.  Proportional  Navigation  Constants  and  Guidance 
Commands 

The  earth  horizontal  and  vertical  line  of  sight 
rates  (o ^  and  which  are  generated  from  the  body  rates 


and  by  the  transformation  are  multiplied  by 

B  B 


a  navigation  constant  and  V^,/g  to  generate  commanded  lateral 
accelerations  in  the  earth  vertical  and  horizontal  planes. 


Target 


A  navigation  constant  of  4.0  in  each  plane  was  found  to 
yield  best  results.  Values  from  3.0  to  4.5  were  evaluated 
with  results  of  that  evaluation  shown  in  Chapter  V.  A  one 
g  bias  is  applied  to  the  vertical  commanded  lateral  accelera¬ 
tion  to  compensate  for  acceleration  due  to  gravity.  Thus 

two  lateral  acceleration  commands  (nv  and  n7  )  are  computed 

c  Zc 

in  the  earth  horizontal  and  vertical  planes. 


n 


Y 


c 


(13) 


n 


Z 


°EL  +  1 


(14) 


c  g 

These  commanded  lateral  accelerations  are  subjected  to  a 
further  transformation  (T^)  in  order  to  generate  commanded 
bank  angle  and  normal  load  factor  which  are  required  by  the 
missile  autopilot. 


C 


tan 


n  =  n,,  cos^+n..  sintfc 

Z  „  a  1 

C  C  C 

Figure  4-7  depicts  the  geometric  relationship  described  by 

the  above  two  equations.  Commanded  bank  angle  (  and 

normal  load  factor  (n  )  are  used  as  inputs  to  the  appropriate 

zc 

autopilot  outer  command  loops  to  provide  closed  loop  bank 
angle  and  load  factor  responses. 


75 


b.  Limits  on  Guidance  Commands 


In  the  attack  phase,  the  following  limits  are 
applied  to  commanded  variables: 

-4 

|  * c  I  —  l  max 

Missile  performance  is  evaluated  for  maximum  roll  rates  of 
50,  100  and  200  degrees  per  second. 

c.  Bank  Angle  Anomoly 

In  the  attack  mode  there  is  no  limit  on 
commanded  bank  angle.  Thus  a  computational  scheme  is 
required  to  insure  that  the  missile  always  rolls  the 
shortest  way  to  a  new  commanded  bank  angle.  This  is  accom¬ 
plished  by  continuously  calculating  the  quantity: 

A  (t)  =  1$  ,(t)  -  -j>(t)  | 

If  A(t)  is  less  than  180  degrees  then  4>c(t)  is  sent  to  the 

bank  angle  autopilot.  If  A(t)  is  greater  than  180  degrees 

then  a  modified  bank  angle  command  ( <fc  )  is  sent  to  the  auto 

M 

pilot  in  accordance  with  the  following  logic: 

<t>  =  <t>„+360°  for  4>  <0 

CM  °  C 

<j>  =  <t>  -360°  for  <t>  >0 

cM  c  C 

d.  Switch  Positions 

The  switches  shown  in  figure  4-6  are  in  the 
following  positions  for  the  attack  phase: 

SW1  Either  position 
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SW2  Position  2 

SW3  Position  1 

SW4  Closed 

e.  Attack  Phase  Summary 

In  the  attack  phase  proportional  navigation  is 
used  in  both  azimuth  and  elevation  to  generate  bank  angle 
and  normal  load  factor  commands  to  the  missile  autopilot. 

The  phase  is  commenced  at  a  range  of  9100  feet  with  missile 
offset  in  heading  from  the  target  at  an  altitude  of  approxi¬ 
mately  250  feet.  Roll  rate  is  limited  to  <i>max*  Bank  angle 
is  not  limited.  Normal  load  factor  is  limited  from  -2  to 
4  g's. 

2.  Ingress  Phase 

The  ingress  phase  commences  with  the  beginning  of 
the  problem  simulation  when  the  missile  is  25,000  feet  from 
the  target  and  initially  heading  north.  The  missile's  Mach 
number  is  0.75  and  the  ship's  velocity  is  21  knots  to  the 
east. 

a.  Altitude  Hold  Methodology  and  Guidance  Commands 

The  missile  is  required  to  maintain  altitude  of 

50  feet  above  sea  level  during  this  phase.  It  accomplishes 

this  by  comparing  present  altitude  with  50  feet  and  multiplying 
K 

the  error  by  ■  to  establish  a  flight  path  angle  command. 

VT 

Kjj  is  assigned  a  value  of  0.3.  This  calculated  flight  path 
angle  command  is  compared  with  the  present  flight  path  angle 
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r 


and  the  difference  or  error  is  multiplied  by  K^and  - 

to  give  an  elevation  line  of  sight  rate,  a  gg.  A  value  of 

0.3  is  assigned  for  Ka.  Thus  n„  is  computed  as  follows: 

Zc 
V, 

(  y  -  y)  +  i 


nzc  =  kh 


_t_ 

g 


The  lateral  acceleration  command  rij  is  computed 


c 

exactly  as  described  in  the  attack  phase  with  bank  angle  and 
roll  rates  limited. 

The  second  transformation  for  the  ingress  phase 
(Tjj)  remains  the  same  form  as  the  attack  phase. 


b. 


applied  to 


n  =  n„  cos  <frf-ny  sin  <J> 
zc  c  ~c 

Limits  on  Guidance  Commands 

In  the  ingress  phase  the  following  limits  are 
commanded  variables: 

-2  <  nz  <4 


|  ^  —  I  4iax  ^ 

The  maximum  bank  angle  allowed  is  60  degrees  and  roll  rate 
is  limited  to  50,100  or  200  degrees  per  second, 
c.  Switch  Positions 

The  switch  positions  as  shown  in  Figure  4-6  are 
as  follows  during  the  ingress  phase: 
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SW1 ,  SW2,  SW3  Position  1 

SW4  Closed 

3.  Offset  Phase 

The  missile  initiates  the  offset  phase  at  17,800  feet 
The  offset  phase  is  subdivided  into  two  sub-phases  --  turn 
and  level. 

a.  Turn  and  Level  Sub-Phases 

The  missile  maintains  an  altitude  of  50  feet 
throughout  the  offset  phase  in  the  same  manner  as  described 
during  the  ingress  phase.  Thus 

nzc  =  Sir  (Y  -Y)  +  1 
s 

The  missile  performs  the  coordinated  turn  by 
using  a  commanded  bank  angle  and  computing  a  normal  load 
factor  as  fellows: 

*  -  60° 

a 

2c  COS  <J) 

At  a  range  of  1 5800  feet  from  the  target  a  turn  of  approxi¬ 
mately  12  degrees  is  accomplished.  The  missile  then  performs 
a  straight  and  level  maneuver  by  commanding  zero  bank  angle 
along  with  the  same  normal  load  factor  as  follows: 

*0  =  0° 


V  ~  COS  <t> 

c 
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b.  Limits  on  Guidance  Commands 


This  commanded  climb  angle  is  compared  with  the  actual  y  and 

vT 

the  error  is  multiplied  by  and  —  to  give  an  elevation 
line  of  sight  rate.  A  value  of  1.0  is  assigned  for  K y 
in  this  simulation.  Thus  n^  is  computed  as  follows: 


<T0li.b  '  Y>  +  1 


The  lateral  acceleration  commanded  n_  is  computed 

^c 

exactly  as  described  in  the  ingress  phase. 

The  second  transformation  T^  is  as  follows: 

$  =  tan"  ^c 

c  - - 

c 


n  =  n7  cosd+nv  sin<j> 

2c  “c  1c 

The  missile  is  free  to  start  rolling  into  the 
target  using  proportional  navigation  during  the  pop  up 
phase  and  thus  aids  in  the  attack  phase  which  follows  the 


pop  up  phase. 

b.  Limits  on  Guidance  Commands 

In  the  pop  up  phase  the  following  limits  are 
applied  to  commanded  variables: 

-2  -  n  -  4 

zc 


<t>  £  <t> 

c  n 

<t>cl  1  I* 


The  maximum  bank  angle  is  limited  to  6C  degrees 
and  the  maximum  roll  rate  is  limited  to  50,  100  or  200  degrees 
per  second. 

c.  Switch  Positions 

The  switch  positions  as  shown  in  Figure  4-6  are 
as  follows  during  the  pop  up  phase: 

SW1  Position  2 

SW2,  SW3  Position  1 

SW4  Closed 

5.  Guidance  Summary  baseline) 

The  following  matrix  indicates  the  horizontal  and 
vertical  controls  during  each  phase  of  missile  flight. 


PHASE 

AZIMUTH  CONTROL 

ELEVATION 

CONTROL 

Ingress 

PN* 

Altitude 

Hold 

Offset 

4>  Command 

Altitude 

Hold 

Popup 

PN 

Flight  Path  Angle  Hold 

Attack 

PN 

PN 

^Proportional  Navigation 

D.  ALTERNATE  GUIDANCE  LAW  DESIGN  SEA-SKIMMER 

The  sea  skimmer  scheme  differs  only  from  the  baseline 
guidance  scheme  in  that  the  offset  and  popup  phases  are 
eliminated.  In  essence  this  scheme  is  composed  of  an  ingress 
phase  and  an  attack  phase. 

1 .  Ingress  Phase 

The  ingress  phase  commences  with  problem  initialization 

which  locates  the  missile  25>000  feet  from  the  target  with 

83. 


identical  initial  conditions  as  the  baseline  guidance  scheme, 
The  vertical  load  factor  is  computed  as 


V„ 


nr 


=  K 


H 


(Yc  -  Y)  +  1 


The  lateral  acceleration  command  is  computed  as 


n. 


_  A 


Vn 


AZ 


g 


AZ 


The  second  transformation  (T^)  is  as  follows: 


n. 


*C  =  tan- 


n„ 


and 


n  =  n„  cos<H-nv  sinb 
Z  L  X 

c  c  c 


Limits  on  guidance  commands  and  switch  positions 
for  the  ingress  phase  are  identical  to  the  ingress  phase 
for  the  baseline  scheme. 

2.  Attack  Phase 

The  attack  phase  commences  at  3000  feet  range  and 
uses  proportional  navigation  in  both  elevation  and  azimuth, 
The  missile  dives  into  the  target  from  50  foot  altitude. 

The  bank  angle  command  and  normal  load  factor  are 
computed  exactly  as  in  the  baseline  scheme  attack  phase 
and  are: 


n. 


=  tan' 


•  1 


nr 
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n  =  n„  cos<t>+nv  sind> 

Z  _  L .  i  _ 

c  c  c 

Limits  on  guidance  commands  and  switch  positions 
are  also  identical  to  the  attack  phase  in  the  baseline 
guidance  scheme. 

3.  Guidance  Summary  (Sea  Skimmer) 

The  following  matrix  shows  the  horizontal  and  vertical 
controls  during  each  phase  of  the  sea  skimmer  guidance 


scheme. 

PHASE 

AZIMUTH  CONTROL 

VERTICAL 

CONTROL 

Ingress 

PN 

Altitude 

Hold 

Attack 

PN 

PN 
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V.  PROBLEM  SIMULATION.  RESULTS,  CONCLUSIONS  AND  RECOMMENDATIONS 


This  chapter  describes  the  computer  simulation  special 
language  used  and  presents  the  results  of  test  runs  in  which 
roll  rate,  ECM  blinking  rate  and  burn  through  range  are  all 
varied  in  both  the  baseline  guidance  scheme  and  the  sea 
skimming  guidance  scheme.  From  these  results  optimum  para¬ 
meters  are  chosen  to  yield  minimum  miss  distance  for  each 
guidance  scheme.  Miss  distance  ranges  are  also  presented 
for  the  baseline  scheme  with  Glint  imposed. 

A.  CSMP  SIMULATION 

The  Continuous  System  Modeling  Program  III  (CSMP  III), 
an  IBM  developed  system,  is  the  modeling  technique  chosen 
for  this  study.  CSMP,  as  discussed  in  Speclchart  and 
Green  {Ref.  4).  is  compatable  with  the  Fortran  IV  language, 
has  six  selectable  integration  methods  (fourth  order  Runge- 
Kutta  is  used  for  this  simulation) ,  and  has  multiple  special 
function  codes  that  can  solve  a  large  variety  of  control 
engineering  problems. 

The  special  functions  used  for  this  simulation  are: 

a.  FUNGEN  -  This  function  extracts  the  value  of  a 
function  of  one  variable  from  a  stored  table  in  memory.  It 
is  used  for  aerodynamic  coefficients. 

b.  TWOVAR  -  This  function  extracts  the  value  of  a 
function  of  two  variables  from  a  stored  table  in  memory.  It 
is  used  for  aerodynamic  coefficients. 


c.  INSW  -  Y  =  INSW  (XI,  X2,  X3 )  assigns  Y  =  X2  for 

X  <0,  Y  =  X^  for  X1  -  0.  It  is  used  to  accomplish  ECM 
blinking  action. 

d.  GAUSS  -  This  function  is  a  random  number  generator 
with  normal  distribution.  It  is  used  to  realize  glint. 

e.  REALPL  -  This  function  simulates  a  first  order 
lag  1/(tS++1)  filter.  It  is  used  in  autopilot  and  seeker 
equations  as  low  pass  filters. 

f.  DEBUG  -  This  is  a  computer  program  aid  that  evaluates 
all  output  parameters  and  variables  at  requested  times. 

In  addition  to  the  special  functions,  trigonometric, 
integral,  derivative  and  square  root  functions  are  used  to 
perform  necessary  steps  with  one  invocation.  CSMP  also  has 
the  ability  to  perform  calculations  in  either  a  sort  section, 
where  the  order  of  computation  is  not  considered,  or  a 
nosort  section,  where  ordered  sequential  fortran  statements 
are  required.  The  program  is  structured  as  shows  below: 

MISSILE  SIMULATION  PROGRAM 

INITIAL 

CONSTANTS 

AERODYNAMIC  COEFFICIENT  TABLES 

DYNAMIC 

N0S0RT 

MISSION  PHASE  LOGIC 
GUIDANCE  EQUATIONS 
ECM  (BLINKING)  EQUATIONS 

SORT 

6  DOF  MISSILE  DYNAMIC  EQUATIONS 
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INNER  LOOP  AUTOPILOTS 
SEEKER  EQUATIONS 

TERMINAL 

OUTPUT  GRAPHS 

The  CSMP  initially  sets  up  the  two  and  three  variable 
coefficients  in  tabular  form  in  memory.  The  run  commences 
with  the  missile  in  steady  state  flight  with  steady  state 
variables  established.  These  steady  state  variables  are 
used  to  enter  the  tables  in  memory  and  extract  coefficients 
which  are  used  to  calculate  aerodynamic  derivatives,  forces, 
moments,  and  other  relations.  From  the  steady  state  flight 
condition  small  perturbations  are  made  and  new  table  entries 
produced.  This  process  is  repeated  for  all  phases  of  flight. 

B.  RESULTS 

Both  the  baseline  guidance  scheme  and  the  sea  skim 
guidance  scheme  were  tested  under  the  same  conditions  and 
quantitative  results  obtained.  All  scenarios  started  from 
the  following  initial  conditions: 

1.  Target  position  -  25>000  feet  north  of  the  missile. 

2.  Target  speed  -  east  at  21  knots. 

3.  Missile  speed  -  north  at  0.75  Mach  number 

4.  Missile  altitude  -  50  feet  above  sea  level. 

The  simulation  was  run  for  31  seconds.  Miss  distance  data 
was  obtained  for  varying  roll  rate  limits,  ECM  blinking 
frequencies,  and  burn  through  ranges.  The  following  matrix 
of  parameters  was  tested: 
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Burn  through  range  (feet) 

2000 
1  500 
1200 
800 
400 

Roll  rate  limit  (degrees  per  second) 

200 
1  00 
50 

ECM  blinking  frequency  (Hertz) 

2 

1 

0.5 

The  resulting  miss  distances  are  shown  in  Tables  5-1 
and  5-2  for  all  runs.  Figure  5-1  through  figure  5-6  are 
the  graphical  results  of  simulation  runs.  Each  run  is 
presented  in  6  consecutive  figures  a  through  f  and  depict 
the  following. 

a.  Commanded  roll  rate  (PCMD)  and  roll  rate  (P)  vs 
time . 

b.  Commanded  bank  angle  (PHICMD)  and  bank  angle  (PHI) 
vs  time. 

c.  Range,  X  Range,  I  Range  and  Z  Range  vs  time. 

d.  Left  and  right  stabilizers  (LSTAB  and  RSTAB)  vs 
time. 

e.  Rudder  position  vs  time. 

f.  Map  of  YMISLE  and  HMISLE  vs  XMISLE. 

The  following  matrix  indicates  the  test  conditions  used 


to  generate  each  figure: 


TABLE  5-1 


Baseline  Guidance  and  Control  Scheme 


Miss 

Distance 

(feet) 

ECM 

Blinking 

Rate 

Roll 

Rate 

Limit 

400 

Burnthrough  Ranges 

800  1200  1500 

.  5  cps 

50 

19.9 

1  5.7 

1  6 

16 

100 

17 

11.3 

6.95 

6.95 

200 

27 

30 

29 

29 

1  cps 

50 

24 

23 

11.3 

9.8 

1  00 

12.6 

13.7 

12.5 

11 

200 

5.4 

3 

4.9 

2.9 

2  cps 

50 

25 

20 

3.9 

6.6 

1  00 

23.7 

1  5.7 

2.2 

1.69 

200 

13.3 

7.9 

1.5 

4.2 

0  cps 

50 

2.5 

100 

4*4 

200 

4.5 

2000 


6.28 

18.16 

26.1 

10.7 

6.9 

6.57 

4.34 

6.8 

6.81 


90 


TABLE  5-2 


Sea-Skimmer  Guidance  and  Control  Scheme 
Miss  Distance  (feet) 


ECM  Roll 
Blinking  Rate 
Rates  Limit 


400 


Burnthrough  Ranges 


800 


1200  1 500 


BASE  GUIDANCE  SCHEME 


5-1 

(a-f) 

NO  ECM,  100°/sec  roll  rate 

5-2 

(a-f) 

ECM  (1Hz),  100°/sec  roll  rate 

5-3 

(a-f) 

ECM  (1Hz),  50°/sec  roll  rate 

5-4 

(a-f) 

ECM  (1Hz).  200°/sec  roll  rate 

SEA  SKIMMER  SCHEME 

5-5 

(a-f) 

NO  ECM,  l00°/sec  roll  rate 

5-6 

(a-f) 

ECM  (1Hz).  100°/'sec  roll  rate 

1 . 

Analysi 

s  of  Baseline  Scheme  Results 

a.  Varying  Burn  Through  Range 

Miss  distances  generally  decrease  as  the  burn 
through  range  increases.  For  runs  beyond  1500  feet  burn 
through  ranges,  the  miss  distances  are  close  enough  to  be 
considered  as  hits  since  most  of  the  miss  distance  is  in  the 
Y  component  of  range  which  has  the  largest  offset  (75  feet) 
applied  during  ECM. 

b.  Varying  ECM  Blinking  Frequency 

The  miss  distance  decreases  as  blinking  rate 
increases  for  the  same  burn  through  range.  Miss  distances 
for  burn  through  ranges  beyond  1200  feet  are  particularly 
good  for  all  three  roll  rates. 

c.  Varying  Roll  Rate  Limit 

Varying  roll  rate  limit  results  in  the  largest 


variation  in  miss  distance.  For  200  degree  per  second 
roll  rate  limit  and  0.5  Hertz  ECM  blinking  frequency  the 
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Figure  5-1 c  Baseline  Response  for  No  ECM,  Roll  Rate  Li 

(RANGE,  XRANGE,  Y RANGE,  ZRANGE  vs  TIME) 
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Baseline  Response  for  ECM  1  CPS,  Roll  Rate  Limit  200  deg/s 
(L  STAB,  R  STAB  vs  TIME) 


( YMISLE,  HMISLE  vs  XMISLE) 


YMISLE,  HMISLE  vs  XMISLE) 


(PCMD,  P  vs  TIME) 


(PHICMD,  PHI  vs  TIME) 


(RANGE,  XRANGE,  Y RANGE ,  ZRANGE  vs  TIME) 
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( YMISLE,  HMISLE  vs  XMISLE) 


miss  distance  remains  high  for  all  burn  through  ranges  out 
to  2000  feet.  But,  for  the  remaining  blinking  frequencies 
and  burn  through  ranges  the  miss  distances  improve  signifi¬ 
cantly  for  increasing  roll  rate  limits.  A  miss  distance  of 
7.9  feet  for  a  burn  through  range  of  800  feet  and  2  Hertz 
blinking  frequency  is  excellent. 

2.  Analysis  of  Sea  Skimmer  Scheme  Results 

The  results  for  the  sea  skimmer  scheme  are  quite 
different  than  the  baseline  scheme. 

Burn  through  range  seems  to  be  the  single  most 
controlling  variable.  Miss  distances  for  a  burn  through 
range  of  1500  feet  are  marginally  satisfactory  for  all  roll 
rate  limits  and  ECM  blinking  rates.  However,  at  a  burn 
through  range  of  2000  feet,  miss  distances  for  varying  roll 
rate  limits  and  ECM  blinking  frequencies  are  extremely  low 
and  considered  quite  good. 

Varying  the  ECM  blinking  rate  produces  little  effect 
on  the  miss  distances.  However,  a  unique  roll  rate  limit 
of  100  degrees  per  second  produces  miss  distances  double  the 
miss  distances  for  roll  rate  limits  of  50  or  200  degrees 
per  second. 

3.  Analysis  of  Glint  Simulation  Results 

Figure  5-7  shows  the  random  introduction  of  glint 
in  the  target  position  and  line  of  sight  rate.  The  miss 
distance  for  the  run  in  this  figure  is  1 5  feet,  but  analysis 
from  other  runs  show  that  glint  increases  the  miss  distance 
30  to  50  percent. 
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4.  Proportional  Navigation  Ratio  Selection 

The  following  matrix  shows  the  result  of  miss  distance 
in  a  no  ECM  environment  for  varying  proportional  navigation 


constant  (PN). 


Miss  Distance  (feet) 

10.6 


Although  a  proportional  navigation  constant  (PN) 
of  4.5  results  in  the  smallest  miss  distance  excessive  over 
shoot  to  step  inputs  results  in  this  proportional  navigation 
constant  as  being  an  unsatisfactory  choice.  A  PN  of  4.0  was 
chosen  based  upon  satisfactory  system  responses  and  the 
snail  miss  distance  achieved. 


G.  CONCLUSIONS  AND  RECOMMENDATIONS 

Based  upon  the  results  of  this  study  it  can  be  concluded 
that  miss  distance  is  a  function  of  roll  rate  limit,  ECM 
blinking  rate,  burn  through  range,  and  the  guidance  scheme 
chosen.  The  guidance  scheme  chosen  is  most  significant 
of  these  factors. 

For  the  baseline  guidance  and  control  scheme,  roll  rate 
limit  is  the  most  significant  variable  affecting  miss 
distance.  One  exception  to  this  is  the  0.5  Hertz  blinking 
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rate  set  of  runs  which  shows  a  miss  distance  improvement 
when  roll  rate  is  increased  from  50  to  100  degrees  per  _second, 
but  yields  poor  miss  distances  for  200  degrees  per  second. 

The  most  significant  variable  affecting  miss  distance 
for  the  sea-skimmer  scheme  is  burn  through  range.  Nevertheless, 
a  1 00-degree-per-second  roll  rate  limit  appears  to  be 
unacceptable  at  all  burn  through  ranges  simulated. 

In  summary,  given  that  the  target  controls  ECM  blinking 
frequency  and  burn  through  range,  the  optimum  selections 
are: 

Baseline  Scheme  -  Roll  rate  limit  =  200  degrees  per  second 

Sea-Skimmer  -  Roll  rate  limit  =  50  or  200  degrees  per 
second 

The  following  additional  studies  are  recommended: 

a.  Further  modeling  with  glint. 

b.  Variation  of  the  direction  and  magnitude  of  the  ECM 
shift. 

c.  Alteration  of  the  decision  variable  used  in  the 
simulation  to  switch  phases  from  range  to  line  of  sight 
rate. 

d.  Variations  in  flight  path  geometry. 

e.  Comparison  of  control  laws  using  skid-to-turn  with 
control  laws  using  bank-to-turn  maneuvering. 
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Side  Slip  Characteristics 
*  vs  3  vs  C 
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Side  Slip  Characteristics 
*  vs  8  vs  C 

nbasic 
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Lift  Characteristics 
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Directional  Control  Characteristics 
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Directional  Control  Characteristics 


*  vs  5  vs  A  C  0;  ) 
nST 
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Pitching  Moment  Characteristics 
*  vs  C, 
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basic 
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Drag  Characteristics 


Drag  Characteris 
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STEADY  STATE  DATA 


2200  LBS  (At  end  of  cruise  flight) 
68.38  Slugs 
12  feet2 
8.485  Feet 

I .  41 4  Feet 
27.8  Slug  feet2 

1507  Slug  feet2 

1  512  Slug  feet2 

II. 7  Slug  feet2 
600  LBS 

0.75 

339  Feet/Sec 


*  PSLVt  =  837  LBS/feet2 

o  o 

- 2 -  =  0.219 

qSSS 

0.0242 


SS 


=  2.65 


SS 


lssk 


cosa 


SS 


=  0.04625  Radian 
243.3  LBS 


3in  °SS 


38.79  Feet/sec 
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838.1  Feet/sec 


USS  =  VTssC0S<tSS 


=  -0.06 


basic 


SS 
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SS 


SS 


=  -1.0  =  0.01 745  Radian 
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APPENDIX  C 


AUTOPILOT  ROOT  LOCUS  PLOTS 


NEW.  AXIS  (UNITS  PEN  INCH)  »  5.0000 
IHOO  AXIS  (UNITS  PEN  INCH)  •  5.0000 

ROOT  LOCUS  OF  PITCH  RATE  LOOP  OF  NORMAL  ACCEL 
COMMAND  FOR  INNER  LOOP  AUJOPILOT  DESIGN 
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4 


f 


NEAL  AXIS  (UNITS  PEA  INCH)  •  5.0000 
I MAC  AXIS (UNITS  PEN  INCH)  •  5.0000 


ROOT  LOCUS  OF  NORMAL  ACCELERATION  COMMANO 
FOR  INNER  LOOP  AUTOPILOT  DESIGN 


r.. 


nem.  axis  (units  nn  inch)  •  1.0000 

I MAG  AXIS (UNITS  fen  INCH)  •  1.0000 


ROOT  LOCUS  OF  NORMAL  ACCELERATION 
COMMANO  FOR  INNER  LOOP  AUTOPILOT  OESIGN 
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ACAL  AXIS  (UNITS  PEA  INCH)  »  5.0000 
IHAC  AXIS (UNITS  PEA  INCH)  «  S.0000 

ROOT  LOCUS  OF  ROLL  RATE  LOOP  OF  BANK  ANGLE 
C0MMANO  FOR  INNER  LOOP  AUTOPILOT  DESIGN 


I MAC  AXIS  (UNITS  PEA  INCH)  •  5.0000 

ROOT  LOCUS  OF  BANK  RNGLECOMMRNO 
FOR  INNER  LOOP  AUTOPILOT  DESIGN 


new.  axis  (units  pen  inchi  «  s.oooo 

I  NRG  AXIS  (UNITS  Ftfl  INCH)  -  S.OOOO 

ROOT  LOCUS  OF  TRW  RATE  LOOP  OF  TURN 
COORDINATOR  FOR  INNER  LOOf  AUTOPILOT  DESIGN 


g 


ACM.  AXIS  (UNITS  PM  INCH)  -  1.0000 
IHAO  AXIS  (UNITS  PM  INCH)  •  1.0000 

ROOT  LOCUS  OF  TRW  RATE  LOOP  OF  TURN 
COORDINATOR  FOR  INNER  LOOP  AUTOPILOT  DESIGN 


59 


I  MAO  AXIS 


ABAC  AXIS  (UNITS  PEA  INCH)  ■  5.0000 
INNS  AXIS (UNITS  PEA  INCH)  •  5.0000 

ROOT  LOCUS  OF  TURN  COORDINATOR 
FOR  INNER  LOOP  AUTOPILOT  DESIGN 


REAL  AXIS  (UNITS  PER  INCH)  «  1.0000 
IHAG  AXIS (UNITS  PER  INCH)  •  1.0000 


ROOT  LOCUS  OF  TURN  COORDINATOR 
FOR  INNER  LOOP  AUTOPILOT  DESIGN 


KM.  RXI3  (UNITS  PER  INCH)  -  5.0000 
1MR0  AXIS (UNITS  PER  INCH)  •  5.0000 


ROOT  LOCUS  OF  FLIGHT  PATH  ANGLE 
COMMAND  FOR  OUTER  LOOP  AUTOPILOT  OESIGN 


NEAL  AXIS  (UNITS  rEft  INCH)  *  1.0000 
I HAG  AXIS (UNITS  PEA  INCH!  -  1.0000 


ROOT  LOCUS  OF  FLIGHT  PATH  ANGLE 
COMMAND  FOR  OUTER  LOOP  AUTOPILOT  DESIGN 


INTO  nm  (UNITS  fEft  INCH!  •  S.0000 

ROOT  LOCUS  OF  ALTITUDE  H0LDCOMMAND 
FOR  -3UTER  LOOP  AUTOPILOT  DESIGN 


KM.  AXIS  (UNITS  PEP  INCH)  ■  1.0000 
I MAC  AXIS (UNITS  PEN  INCH)  -  1.0000 


ROOT  LOCUS  OF  ALTITUDE  HOLDCOMMRND 
FOR  OUTER  LOOP  AUTOPILOT  DESIGN 


-10. 00 


1.00 


■0.00 


00 


-o.oa 


-i.oo  -if.  00 
REAL  RXIS 


-3.00 


-l 
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APPENDIX  D 


BASELINE  GUIDANCE  LAW  SIMULATION 
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APPENDIX  E 


SEA  SKIMMER  GUIDANCE  LAW  SIMULATION 
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